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FOREWORD 


The analysis in this publication was intended to be an analytical tool to 
assist a study of latchups in support of the Magellan space project. The study 
is not yet completed but it is expected that the analytical methods that have 
been developed so far may have other useful applications in the analysis of 
integrated circuits and therefore they are published here as a separate report. 
The analysis applies to the steady state current distribution in a three 
dimensional integrated circuit. A device physics approach — based on a 
perturbation method — rather than an equivalent lumped circuit approach is 
used. The perturbation method allows the various currents to be expressed in 
terms of "elementary solutions" which are solutions to very simple boundary 
value problems. Work is continuing and it is hoped that there will be a Part 2 
report which applies the methods to the steady state characteristics of latched 
states. Examples which illustrate the use of the analytical methods are not 
given in this Part 1 report since this is expected to be discussed in detail in 
Part 2. Comparisons between theoretical predictions and measurements are also 
deferred to Part 2. The most obvious limitation of the present version of the 
analysis is that all depletion region boundary surfaces are treated as 
equipotential surfaces and this is why the report is subtitled, "A Simple 
Steady State Theory." This may be an adequate approximation in some 
applications but it is an obvious weakness in the theory when applied to 
latched states. This issue is expected to be confronted in Part 2. 
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ABSTRACT 


The analysis applies to the steady state current distribution in a three 
dimensional integrated circuit. A device physics approach — based on a 
perturbation method — rather than an equivalent lumped circuit approach is 
used. The perturbation method allows the various currents to be expressed in 
terms of "elementary solutions" which are solutions to very simple boundary 
value problems. 

The most obvious limitation of the present version of the analysis is 
that all depletion region boundary surfaces are treated as equipotential 
surfaces and this is why the report is subtitled, "A Simple Steady State 
Theory." This may be an adequate approximation in some applications but it is 
an obvious weakness in the theory when applied to latched states. 

Examples that illustrate the use of these analytical methods are not 
given in this report because they will be presented in detail in a future 
report. 
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SECTION 1 


INTRODUCTION 


The analysis in this report was intended to be an analytical tool to 
assist a study of latchups in support of the Magellan space project. An 
investigation of the application of this method of analysis to the steady 
state characteristics of latched states has not yet been completed. A time 
dependent version of this analytical method has also not yet been completed 
but work is continuing in these areas. So far, the analytical method, 
applicable to the steady state current distribution in a three dimensional 
integrated circuit, is reasonably well developed. Some work remains, as will 
be pointed out below, before it will be effective in the study of the steady 
state characteristics of latched states but it is expected that this analytical 
method will have many useful applications in its present stage of development. 
Therefore, it has been published here as a separate report. It is anticipated 
that there will be a Part 2 report which discusses the application of this 
theory to the steady state characteristics of latched states. Examples which 
demonstrate the use of the analytical method are not included in Part 1 since 
this will be treated in detail in Part 2. Comparisons between theoretical 
predictions and measurement are also deferred to Part 2. The most obvious 
limitation of the present version of the analysis is that all depletion region 
boundary surfaces are treated as equipotential surfaces and this is why this 
report is titled "...Part 1: A Simple Steady State Theory." This 

approximation is probably adequate in many applications but it is a definite 
weakness in the theory when applied to latched states. If there is not a more 
elegant way to deal with the boundary values at a depletion region, we could at 
least partition it into sections with each section treated as an equipotential 
surface. This would be an improvement over treating the entire depletion 
region surface as an equipotential. However, such a partitioning would require 
that we correctly represent the tangential currents that flow inside the 
depletion region between adjacent sections and it is not obvious at this time 
how to do this. This issue will be confronted in Part 2. Finally, it is 
hoped that there will eventually be a Part 3 report which discusses transient 
conditions. 


1-1 



The customary method of analyzing the behavior of a real (three 
dimensional) integrated circuit is to simulate the device with an equivalent 
lumped element circuit. A problem with using a lumped circuit to simulate a 
distributed system is that the interactions between various regions in a device 
are not the same as they would be if the various regions were separated and 
connected by wires. It is sometimes difficult to determine representative 
values for the circuit elements used in an equivalent circuit. This is 
especially true for the gains of parasitic transistors when the device is in a 
latched state [1]. These problems can be avoided by treating the device as a 
distributed system, which is what it is. A difficulty with this approach is 
that the governing equations are nonlinear and difficult to solve unless some 
approximations can be made. An approximation proposed in this report is to 
treat the electric field as a perturbation, that can be treated by first order 
perturbation theory, in the quasineutral regions of the device. This produces 
a system of linear differential equations in the quasineutral regions and a 
superposition principle will be used to construct solutions out of "elementary” 
solutions. The elementary solutions satisfy very simple boundary value 
problems. Boundary conditions, which represent the effects of various kinds 
of junctions (e.g., depletion regions, high-low junctions, etc.), will have to 
be satisfied and this will result in a system of nonlinear equations. However, 
these equations will be algebraic rather than differential equations. The 
elementary solutions will be used to construct various constants (the constants 
will be surface or volume integrals of these solutions). These constants will 
be analogous to device parameters (e.g., resistances) used in an equivalent 
circuit simulation and they can be evaluated either experimentally or by 
analytically or numerically solving very simple boundary value problems. With 
the constants regarded as known quantities, the algebraic equations which 
represent the boundary conditions can be solved. 
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SECTION 2 


APPROXIMATION FOR THE ELECTRIC FIELD 


Consider a quasineutral region which is denoted by R and assumed to be 
homogeneous. The well-known equations for the electron and hole currents are 
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eD [- + p y ] 

r v T 


(la) 
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eD 
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+ n 
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where V T = KT/e is the thermal voltage and the other symbols have the 
obvious interpretations. Since R is a quasineutral region, we assume that the 
"applied" electric fields are much stronger than the fields associated with 
space charge separations in R and we make the approximation that the field ? 
in (1) is the applied field. This applied field is associated with surface 
charges induced by an external power supply rather than charges on the 
interior of the homogeneous region R. So the first approximation that is made 
is the assumption that the electric field in (1) has a zero divergence in R. 
This approximation (which obviously does not apply inside transition regions) 
is nothing new. Current induced voltage drops are typically solved by 
assuming ohmic drops and solving Laplace's equation [2]. This is equivalent 
to the assumption of zero charge densities in the quasineutral regions. 
However, "applied" field has not yet been precisely defined. The only 
quantitative statement made so far is that it has a zero divergence in R. It 
is not completely specified until boundary values are specified for the 
potential that the field is derived from. Approximations that will be used 
for the boundary values will be discussed in the next section but, ideally, 
the potential that the electric field is derived from is to be the same as the 
"applied potential" which is defined more rigorously in discussions of 
quasi-Fermi levels. Working with this kind of potential makes it unnecessary 
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to subtract out the built-in equilibrium potentials associated with each 
transition region in order to compute terminal voltages. 

An additional approximation will be the use of first order perturbation 
theory, with the electric field regarded as the perturbation. There are 
intuitive reasons why we might expect first order perturbation theory to work, 
but the reasons applicable to majority carriers are different than the reasons 
applicable to minority carriers. In the case of minority carriers, the drift 
current is typically small compared to the diffusion current since most of the 
drift current is carried by the majority carriers. In many situations, the 
drift term in the minority current can be completely neglected [3]. In 
situations where it can't be completely neglected, it is reasonable to expect 
that it can be treated as a small perturbation. In the case of majority 
carriers, the drift term can't be treated as small since it is often the 
dominant term. However, the success of first order perturbation theory does 
not necessarily require that the perturbation be small, it does require that 
the higher order terms in the expansion of the unknown function be small 
compared to the linear terms. The customary use of simple resistances to 
calculate current induced voltage drops suggests that the majority current is 
approximately a linear function of the applied electric field (this is not the 
same as saying that currents are linear functions of terminal voltages. 
Nonlinearities result from boundary values at various junctions) and this 
suggests that first order expansions may be adequate. 

Although the electric field in (1) will be approximated as being the 
applied field, the field is still an unknown (if it were known, perturbation 
methods would be unnecessary because the equations would already be linear) 
because some of the boundary values that control it are unknown. In fact, 
these boundary values will be the parameters used in the perturbation 
expansion. Therefore, boundary values should be discussed before the details 
of the perturbation method are discussed. 
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SECTION 3 


BOUNDARY CONDITIONS 


The boundary surface of the region R is assumed to be composed of discrete 
sections where any given section is either an insulated surface (oxide layer 
or semiconductor to free space junction), depletion region surface, metal to 
semiconductor contact, or a high-low junction. 

For each boundary surface that separates two interacting regions, the 
boundary conditions will be characterized in two different ways. The first 
characterization states the general class (e.g., Dirichlet, Neumann or 
mixed) that the boundary conditions fall in. A typical example of such a 
characterization is the statement that the potential and carrier densities are 
assumed constant on a given surface. The values of these constants, i.e., the 
boundary values, are often unknowns and must be adjusted to satisfy the second 
characterization (discussed below). The first characterization will be used 
when working with the differential equations to express solutions in terms of 
the unknown boundary values. Later, after the equations have been reduced to 
algebraic equations, the second characterization is to be used to complete the 
solution. The second characterization relates boundary values on opposite 
sides of a given surface. If all potentials were known quantities and we are 
interested in solving for majority and minority currents and/or densities, 
four equations are needed for each boundary surface. In their most simple 
form, the equations would be specifications of both kinds of carrier densities 
on both sides of the surface (in this context, "surface" refers to a junction 
even though it has a finite thickness so carrier densities can be different on 
different sides of the surface). More generally, the equations could include 
a relationship between carrier densities on opposite sides of the surface and 
a relationship between current densities on opposite sides of the surface. 

Since the potentials are also unknowns that must be solved, the boundary 
conditions must also contain enough information to specify the potential on 
both sides of the boundary surface. Altogether we will need six equations for 
each boundary surface that separates interacting regions. The surfaces of 
this type that are discussed in this report are high-low junctions and 
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depletion regions. Surfaces that are not of this type are insulated surfaces 
and metal to semiconductor contacts. For these surfaces it is only necessary 
to solve the equations on one side of the surface. Therefore, only three 
boundary conditions (applicable to the relevant side of the surface) must be 
given. This information is built into the first characterization of the 
boundary conditions so only the first characterization is needed for such 
surfaces . 

A. INSULATED SURFACES 

To find the boundary conditions for the potential, we use the same kind 
of reasoning that would be used if Ohm's law could be applied to the region R. 
There will be surface charges on the boundary surfaces of R associated with 
the applied voltages. If we look within the surface charge layer we would see 
a normal component of the gradient of the charge density and we would see drift 
currents and diffusion currents opposing each other and Ohm's law would not 
apply for any kind of material. Furthermore, the potential does not satisfy 
Laplace's equation here. But if we look just below this surface charge layer 
(in this context, "below" means locations on the interior of R that are very 
close to the surface charge layer) and assume Ohm's law to apply, we would 
conclude that the electric field has a zero normal component at these 
locations. Therefore, the boundary condition that will be used for the 
applied potential is 

fy(X) • K a 0 (2) 

A 

where K is the unit vector perpendicular to the surface. Note that the 
boundary condition (2) is expected to provide an approximation for the 
potential on the interior of R but this region excludes the surface charge 
layers. The currents calculated from this potential exclude the tangential 
currents flowing through the surface charge layers. In particular, the 
current under the insulated gate of a MOSFET is excluded. Since this current 
can be important, it must be included as a separate term. 
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Boundary conditions for the carrier densities will make use of the 
familiar recombination velocity which gives 


— ^ • K = S(p - p )=- — • K for n type 
e * r o e c 


J J 

^ • K=S(n-n Q ) = - ^ • K for p type 


where n and p are the equilibrium carrier densities in R and K is the 
o r o 

outer normal unit vector and S is the recombination velocity. The first 
equations are used if R is composed of n type material and the second are used 
if R is p type. Since the drift currents are assumed to be tangential to the 
boundary surface, the above equations become 


-0^ ^ p • K = S(p - p Q ) = — ^ n • K for n type 
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-D ^p»K=S(n-n)=-D ^ n • K for p type. 

p r On tr Jtr 
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B. DEPLETION REGIONS 


The two boundary surfaces of a depletion region are each assumed to be 
a constant potential, constant n, constant p surface. Using these 
approximations, boundary conditions are of the form 


4>(x) 

= <t> 

= constant 

n(X) 

= n 

= constant 

p(X) 

= P 

= constant 


(4a) 

(4b) 

(4c) 


The boundary values are not independent and at some stage in the analysis they 
must be adjusted to satisfy certain conditions. If the junction is abrupt and 
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if recombination in the space charge region can be neglected, it is shown in 
Appendix A that the conditions are 
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where 


fl = exp [(^ - <|> )/V T ] 


( 6 ) 
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N and N. are the donor density on the n side and the acceptor density on 
Dn Ap 

the p side. The hole and electron carrier densities on the p and n side are 

denoted by p p , n and n while the equilibrium values of these densities 
p n p n 

n and n . The potentials on the two sides are and (j> n . 
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The equilibrium densities can be solved from 


n. 


on 


= if- ex P t^/v T ] 


(7a) 


Dn 


3-4 



(7b) 
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where is the intrinsic carrier density (the density in a pure material) 

and AV and AV are the values of AV in the n material and p material 
n p go 

respectively with AV^ q an empirically determined function of the doping 
concentration. Empirical information on AV^ q has been provided by Slotboom 
and De Graaff [4]. 


To the equations we must add 
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C. HIGH-LOW JUNCTIONS 

The two boundary surfaces of a high-low junction are each assumed to be a 
constant potential, constant n, constant p surface so the boundary conditions 
are of the form (4). 

As with the depletion region, the boundary values must satisfy certain 
conditions. If the junction is abrupt and if recombination in the space 
charge region can be neglected, it is shown in Appendix A that these 
conditions for an n-n + junction are 
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The symbols are the same as with the p n junction except that subscripts k and 

h, which refer to the lightly doped and heavily doped regions, replace p and n. 

AV refers only to the heavily doped side. For a p-p + junction, the equations 
go 

are 
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In both cases we have 
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D. METAL TO SEMICONDUCTOR CONTACT (OHMIC) 

Metal to semiconductor contacts will be assumed to be equipotential 
surfaces with an infinite recombination velocity. Boundary conditions are 
of the form 
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If the contact resistance can be neglected, <J>' in (16a) is the potential of 
the contact lead. If this resistance must be included, it can be simulated by 
placing an imaginary resistor between the contact and the contact lead. In 
this case, <J>' would be the potential on the contact side of the resistor. 
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SECTION 4 


CALCULATION OF CARRIER DENSITIES IN TERMS OF BOUNDARY VALUES 


This section will relate carrier densities to the boundary values by 
assuming boundary conditions to be of the form (2), (3) or (4). Some boundary 
values will usually be unknown and must be adjusted to satisfy equations like 
(5), (8), (9) and (15) but these adjustments are to be made after algebraic 
equations that relate currents to boundary values are obtained. 

To be definite, suppose the region R is an n type material. Note that 
all boundary surfaces except the insulated surfaces specify Dirichlet boundary 
conditions for <j>, p and n. Imagine the boundary surface of R to be the union 
of surface sections with the i ^ section denoted by A^. Let INS and OTH be 
index sets such that: 

If ielNS then A^ is an insulated surface with boundary conditions 
$ <j> • K = 0, -D p $p»K = S £ (p - p Q ), -D n $ n • K = S^p - p Q ). 

If ieOTH then is any surface other than an insulated surface and 
has boundary conditions 4> = \L, n = n^, p = p^. 

In this section, n^ and p^ refer to boundary values on the i th surface rather 
than intrinsic densities. 

Define 


N = n- n, P = p - p . 

o’ r r o 

Using (1) and the equations [3] 
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gives 


(17a) 
(17b) 

1/2 1/2 

where L = (D x) , L = (D x) are the diffusion lengths for holes and 
p p n n 

electrons and x is the recombination time. The boundary conditions can be 
expressed as 


V 2 P - I • $P = \ 
T L 

P 

v 2 n + J-I«^n = ^ 

V T L 

n 


If i e INS and X on then 

$ <J> • K 
$ P • K 

$ N • K 

If i e OTH and X on A i then 

<l> » v i 

P * P. 

1 

N = N. . 

i 



For each j e OTH, define a function <Jk(X) by the conditions 
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J 


(18a) 

(18b) 

(18c) 


(19a) 

(19b) 

(19c) 


(20a) 


If i e INS and X on A £ then $ $ • K = 0 
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= 5. . 
ij 


(20c) 


If ieOTH and X on A. then d> . 

i J 


where S.. is the Kronecker delta. Note that the sum 
ij 


I V <m£) 

jeOTH J J 


satisfies Laplace's equation in R and the boundary conditions imposed on the 
potential 4>. It is well known that solutions to such problems are unique so 
we have 


4>(x) = l v <t> (x) . 

jeOTH J J 


( 21 ) 


The (f^’s are not all independent. The sum 

l <MX) 

jeOTH J 


is one on every noninsulated surface and has a zero normal gradient on the 
insulated surfaces. The constant value 1 satisfies Laplace's equation and also 
has these properties which implies 


1 <MX) = 1 . (22) 

jeOTH 3 

Constraint (22) is consistent with the well known fact that any potential can 
be chosen as a reference and set equal to zero, which eliminates one of the 
functions that must be included on the right side of (21). If none of the 
potentials are set equal to zero, (22) can still be used to reduce the number 
of boundary value problems that must be solved in order to obtain all of the 
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Using (21), the electric field is given by 


E = <t> = - l V. ^ 
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so (17) becomes 
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Equations (24) and the boundary conditions (18) and (19) define P and N as a 
function of the V^'s, i.e.. 


P = P(X, V^, N = N(X, Vp . 


It is assumed that the functions can be expressed in a series of the form 

P(X, V.) = p (o) (X) + l p (X) V. + 0 (2) (25a) 

1 jeOTH J J 1 

N(X, V.) = n (o) (X) + l n.(X) V + 0,(2) (25b) 

1 jeOTH J J 

where 0^(2) and 0^(2) represent terms of second or higher order in the V^'s. 
Parentheses are used to emphasize that the superscripts are indices rather 
than exponents. Substituting (25) into (24) and equating coefficients to like 
powers of the V^’s yield 
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By using the identity 


7 2 (fg) = 2$f • $g + f7 2 g + g7 2 f 


together with 7 <f> = 0 and (26), we can express (27) as 
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To obtain boundary conditions satisfied by the n's and p's, first consider 
an insulated surface A^. From (18b) and (25a) we have, for X on A^, 
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This is to be satisfied identically in the V's which gives 






P 


with analogous equations for the n's. Note that (25a) and the requirement 
that (19b) be satisfied identically in the V^'s give 


= P. , p . =0 
r i J 


on the noninsulated boundary surfaces. The boundary conditions are summarized 
below. 
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(33b) 
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For each jeOTH, define the function ^ by the conditions: 
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If icOTH and & on A. then 
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It is evident that the sum 


2 P. 

jeOTH J J 


satisfies the same field equation and boundary conditions as It is shown 

in Appendix B that solutions to such problems are unique which implies 
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From (30) it is seen that 


P j + 2 V^ P <0> 
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satisfies the same field equation as y . . Note 


? + 2^ p<0) *i ] * * - s • * + Wj p<0) ?p<0> • i ■ 


Using (32b), (20b) and (32a) shows that on the insulated surfaces. 


$ tPj + P <0> *jl - K - 


-S.x 

1 


Pa ~ 


d> . S.x p 

T„ X r 


(0) 


2V, 


T L 


S.x . / \ 

[p j + 2V^ P V 


which is the same boundary condition satisfied by Yj on the insulated surfaces. 
Also, on a noninsulated surface A^, we have from (33a), (33b) and (20c) that 


1 (o) , i f 1 r 

P -4 + 947 P *4 - 947 “i 4 - 947 1 -4 


p. 


j 2v_ 


r j 2 V t ij 2 V t ij 


which are the same boundary conditions satisfied by 


P. 


2V t f j 


We therefore have 


1 (o) , 

P 4 + ^TP *A = 


P. 


j 2V_ 


j 2 V t T j 


or, using (37), 


’i = 2V~ [P i *1 " +1 X P K V 
J 2V T j J J Ke0TH K K 


(38) 


4-8 



Now consider the expression 


(o) 


L 2 , . 

-f P (0) . 

L 

n 


It is seen from (31) that this expression satisfies Laplace's equation inside 
R and from (32a) and (32c) it is seen that on the insulated surfaces. 


? tn <0) - P (o) ] • K = 0 


n 


From (33a) and (33c) it is seen that on a noninsulated surface we have 


(o) 


L 2 

_e 


(o) 


= N. - — r P. 
i .2 i 


The expression 


l (N. - 
jeOTH J 


j* V 

n 


also has these properties so we conclude 


» 


L 2 , , 

p (o) 

L 2 

n 


= 1 (N • - 

jeOTH J 


n 
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or, using (37), 


( \ 1,2 

■ <o) - i l 

IT jeOTH J J 


+ l (N. 
jeOTH J 


i V *j 


- 2 

jeOTH 


[N. <j>. 
3 J 


n 


♦j)] • 


( 39 ) 


Equations (37), (38) and (39) solve for p^°\ and t/°^ in terms of the 
elementary functions and <J k. If the second term on the left side of (29) 
contained a plus sign instead of a minus sign, rij would be easy to solve using 
the same basic methods that were used to solve for the other quantities. As 
it is, n. is not easy to solve. However, it is unlikely that we would be 
interested in the majority carrier density away from the boundary surfaces 
except for the purpose of calculating the gradient in order to obtain the 
diffusion current. Fortunately, as will be shown in the next section, it is 
not necessary to solve (29) to do this. 

The important equations in this section and the analogous equations that 
apply to a p type material are listed in the summary table which follows the 
next section. 
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SECTION 5 


EXPRESSING CURRENTS IN TERMS OF THE BOUNDARY VALUES 


The previous section related carrier densities to the boundary values 
(some of which will usually be unknown) and the elementary functions <j>., y.. 
This section will do a similar thing for the currents through a given 
noninsulated surface. The minority carrier recombination current at an 
insulated surface is also easy to express in terms of the elementary functions 
(this also solves for the majority carrier recombination current since the net 
current is zero) but these currents have no immedidate use in the analysis 
presented here. Therefore, we will confine our attention to the noninsulated 
boundary surfaces. 

For notational brevity in the equations that will follow, some constants 

will now be defined. For each i and i in OTH, define Q. . and R. . by 

1 3 ij 





ds 


(40) 


R. . = /. $ <t>. • ds 

ij 


(41) 


where the integrals are surface integrals on the surface and the unit 
normal vector is in the direction of K, i.e., it is directed outwards from the 
region R. 


R. . is easily recognized as being proportional to the current in an Ohm's 
law medium, through the surface A^, divided by the potential on the surface A^ 
when all noninsulated surfaces except A^ are grounded. R„ is therefore 
inversely proportional to an electrical resistance when the region R contains 
an ohmic material. Q„ is proportional to the current in a pure diffusive 
medium (a medium that supports diffusion currents but not drift currents) 
through the surface A^ divided by the excess carrier density on the surface 
Aj with the excess carrier density zero on all other noninsulated surfaces. 
Therefore, can be thought of as the reciprocal of a diffusion resistance. 


5-1 



The constants Q.^ and R„ 
constraints connecting them. 


are not all independent , there are several 
For example, using (36) gives 


Q. 4 ~ Ja ^ ^ # d ® = l ^ • ds =» f ^ tj; • ds 

lj A i J keOTH \ 1 J 1 3 


- I 

kelNS 


V* 



• ds . 


But (35) gives 


/ • dl . 


L“ \ 
P 


ds 


and Green's first identity gives 


/ $ <|i. • ds = /„ [»|». V 2 . + $ y. • $ * ds 

1 J K 1 J 1 J 


Using (34) gives 


||f # ||f # 

$ t £ $ * ds = X R [-iy- 1 + ^ t £ • ^ tj] d a 


so that 


Qij = S R [_ V^ + * * * *J d J 

L P 


V 


• 7 tji ] d X + 2 T J*A. ♦i ♦i 

1 J kelNS IT Tc 3 


ds 


The expression on the right of the above equation is symmetric in i, j which 
proves the reciprocity relation 



( 42 ) 
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A simpler version of the above proof will prove the well known reciprocity 
relation for the electrostatic potential which, in this notation, is written as 




(43) 


Another constraint can be obtained by taking the surface integral of the 
gradient of both sides of (22) on any noninsulated surface A^. The result is 


I J. $ <j>. • dt = 0 

jeOTH i J 


or 


l R. . = 0 
jeOTH 1J 


for any ieOTH 


(44) 


Equation (44) could also have been deduced by using (20b) and recognizing that 
has a zero surface integral on a closed surface. 


The hole current directed outward through a noninsulated surface is 


given by 


+ + * ( p + p Q ) * + 

I = S. J • ds = e D /. [-V P + — E] • ds 

P £ A. p p A £ V t 


Using (19b) gives 


p i * -► p i + p 0 ± + 

5T - -'a. * p • d * ♦ j 'a. J * d6 

pi T i 


To first order in the V's we set 


p = p ( o) + l V, p. . 


jeOTH 


3 J 
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Using (23), (37) and (38) gives 


S A. ?P 


* 

ds 


V. 


= 1 p ; J'a $ * ds [1 + 2V~] - l l 


L . j A. 
J i 


V. p i, 

J is 

2V 

j k T 


; A. *le ? ' dl 

1 J 


and 


- I 1 !±_!k x * ^ u, . dl 

j k 2 V t J A. 9 j V k 


/ I • dl = - l V / $ <>, • ds . 

A i j J i J 


Using (20c), (36), (AO) and (Al), the equations reduce to 


X A. ^ P * dS " ? P j C1 + j 2V T ] Q ij P i ? 2V^ R ij 


(A5 ) 


'a.* 

x 


ds 


-I V. R. . 
j J iJ 


(A6) 


so that 


ar - - p j [1 + 

P J 


V. - V. P. V. 

J 2V t 1] Q ij " (p o + F' ) ^ R ij' 


(A7) 


It is to be understood that the summations are over noninsulated surfaces so 
the subscripts OTH to the summation symbols were left out. 


To express the majority current (electrons in this case) in terms of the 
boundary values, it is helpful to define the following additional geometric 
constants . 


5 -A 



For i, j, k e OTH, define 


S ijk - /r *i * tj • ? * k 0 3 X 


( 48 ) 


T ijk * X R *i ? ? Ok d3x • 


(49) 


A relationship between T. and R. . can be obtained as follows. Green's 

ijk ij 

theorem can be used to obtain 


X R [4> k V 2 (<j>. <fr.) - (*. $.) V 2 $ fc ] d 3 X 


- X [<|> k $ (<|> i 4>j) - (<t>. 4>j) $ 4> k ] • ds 


or 


2X r 4* k ^ ^ ♦j d 3 x = # [ 4 > k <j). ^ 4> j + <j> k 4>j $ 4>i - 4^ 4>j $ 4> k J • ds 

= 5 ik X* ^ 4>. • ds + X A $ <t>,. • ds 


\ 


- S. . J\ $ * ds • 

ii A. T k 
J 1 


From (41) and (49) we have 


2T, . • =5 5,. R M + 6 M R m — 6. . R m 
kij ik jk jk lk ij lk 


or 


T ijk = 2 tS ij R jk + 6 ik R ij “ 8 jk R ik ] * 


(50) 
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One constraint between the constants S.^ k can be obtained by using the 
divergence theorem with 


§ 4k • y. $ 4>v • ds = S i; j { ^ 4> k * ds = S. . R. 


r i T j ' V. 


A. 

1 


ij ik 


to get 

X r. $ 4>i • $ 4> k d 3 x + X 4». $ r. • $ 4> k d 3 x = R. k 


Interchanging i and k and subtracting equations and then using (48) yield 


S kji “ S ijk + lS jk 6 ij ] R ik * 


(51) 


Another constraint comes from (22) and (48) which immediately give 


i = 0 • 


(52) 


Equations (51) and (52) are the only independent constraints that this author 
has been able to find. If the reader is aware of any other constraints that 
are independent of (51) and (52), this author would appreciate being informed. 


The electron current through the surface is given by 


I = f. j • ds 
n. A. J n 

i i 


or 


n. 


^ = X A $ N • ds + 7T- X A (N + n ) f • ds = X A ^ N 
eD A. V. A, o A. 

n l T l l 


• ds 


(N i *■.> * 

* ~ X A. E ' ds • 
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Using (23) and (41) gives 


5T - V ? N * dt — (N i * n o ) ? ^ R ij • 


(53) 


j T 


Using (25b), (37), (40) and (41) gives the first order equality 


S A $ N • ds 
A. 

l 


f. $ • ds + v. J\ $ n* • ds 

•i j 2 A i 2 


L 2 

r-E 


j p i Q i: * (n j 


-? P.) R. .] + J V. J\ ^ n. • ds . 
T 2 J i 3 L . J A. 'j 

n 2 1 


(54) 


Note 


X A $ Hj * ds = § <j> i $ Hj • ds - X INS 4» i ^ Hj * ds 


where the subscript INS means to integrate on the insulated surfaces. But 
(32b) and (32d) show that on the insulated surfaces. 


* n. • k 

3 


i ’ 'j • ; 


n 


so that 


J. ^ n . • ds 

i 2 


/ ^ Hj • ds 


L 2 

^2 J'lNS *i ? P j 

n 


ds 


Since rij ^ is zero on all surfaces, the equation can be written as 
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X A ^ n . • ds = X[<(>. ^ n. - n. ^ 4> - ] 
j i j j i 


^ l 2 + 

d ® ■ 72 J'lNS *i ^ Pj • d * 

L 

n 


Similar steps produce 


^A. ^ Pj • d ® " ^ Pj - Pj ^ 4> i ] • ds - X INS <l> i $ Pj • ds 


Combining the above equations yields 


X ^ H • • ds 
A i 3 


L 2 

-f X A ^ p. • ds = #[♦. ^ n. - 

L 2 A i J 1 3 

n 


n . ^ <t> • 1 • ds 
J i 


~2 ^ Pa 

L J 

n 


p . ^ <b . ] • ds 
J T i 


Green's theorem together with (27) and (29) give 


2 2 

'a. ? • dJ - 72 J'a. ? Pj • • 'a *i 1,2 "j - 72 72 fjl <j3x 

1 Ll L 

n n 


= V - ^ ^ + ~2 ^ P^ * ^ d3x 

T L 


Using (37), (39), (48) and (49), the equation becomes 


X A $ Ha * ds 

A i J 


L 2 

" “2 X A ^ p j * d ® = v" ~2 ^ P k S ikj 

L 2 1 3 V T L 2 k ■ 1KJ 

n n 


*h ( \ 

Ik 


_£ p ) T 

.2 V ijk 


n 
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But, from (38), 


S k. ? p j • d * * 2Vl J 'A. [P j 5 *j - ?< *j l P A>1 ' dt 

1 X 1 1C 

= rj- {P. Q. . - P. R. . - 8. . I P. Q..} 
2V„ ' j x ij l ij ij f k x ik' 


so the equation becomes 


X A. ? n j 

1 j 


ds 


Lj { ip Q 
V T l 2 *2 y ij 
n 


\ P. R. . 

2 l ij 


2 5 ij X P lc ^ik 


+ i V- 

p 


- x p k T i:k * 2 X p k W 


Combining this with (53) and (54) and using (50) finally produce 


n. 
i 

eD 

n 


N. D V . D V. - V. 

-<F + % + D 2 P i> ? “ij + ? (N j - D 2 V [1 + _J 2 vT J 1 R 
n jT J j J n J T 


ij 


D V. - V. D V. 

+ D 2 ? P j [1 + J 2V ] Q ij + 2 D 2 J P k S ikj 
nj T J n j ,k T J 


(55) 


Equation (55) is a cumbersome way to express the electron current. A 
simpler expression can be obtained by making additional approximations which 
are based on the excess electron density being very small compared to the 
equilibrium density. In order to make use of this, it is necessary to be able 
to identify the quantities that are very small. To first order, the drift 
current is given by 


eD eL 

(J )drif t = (N + n ) E = (N + n ) E 
n V T o xV T o 


eL 2 eL 2 (N. + n ) 

(I )drift = J\ (N + n ) E • ds = — — — ^ — f . E • ds . 

n i xV T A i o xV T A. 
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It is obvious that 


eL 2 N. 
n 


tV„ 


1 s 

J A 


E 


ds 


-eL 


tV, 


- ». I V, J. $ ♦, 


A. 

1 


ds 


is a very small contribution to the drift current. It is reasonable to expect 
that this quantity is negligible compared to the total electron current I n ^. 

It is also reasonable to expect that if we sum the terms over i, the result 


-eL 


] v j x s \ ? *j • dl 


is small compared to I . Note that 

i 


l N, s $ (|). • ds = $ n (o) $ <|>. • ds = /_ ^ n (o) • ^ <|>. d 3 x 

k k \ 2 J R J 


(56) 


so the expression (56) becomes 



} V. 




which can also be identified as a small quantity. This should be small 

compared to I regardless of how the V-j's are chosen so the terms 
n • j 

l 



v. /_ $ n (o) • $ d 3 x 

j R j 


(57) 


should individually be small. Note that since <J>^ satisfies Laplace's equation 
in R, its maximum and minimum values are on the boundary surface of R. The 
boundary values are zero and one so (|>^ should be between zero and one 
everywhere in R (this isn't strictly correct because the boundary values were 
not specified on the insulated surfaces, but it is still reasonable to expect 
<{>k to be on the order of one or less everywhere). But if the expression (57) 
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is small and if 4^ is on the order of one or less, it is reasonable to expect 
that the quantity 


^ V ♦„ « n <0) • ? d 3 x 

is also small compared to I . Using (39), (48) and (49), the above 

i 

expression can be written as 


eL 2 V. 
n 3 

tV m 


ft l \ S kHj * £ (N k - V W 

n n 


But if the above expression is small for independently chosen boundary values, 
the individual terms must be small, i.e.. 



V. 

T 


P a S kHj 


and 


(58) 


eL 2 V. 
n 1 

TV m 


(N« - 


L 2 

~2 V 

L 

n 


kfi-j 


(59) 


can each be identified as small quantities. Expression (56) and equation (41) 
produce other small quantities which are 


eL 

-=r N * v * R - • 

tV t i j ij 


(60) 


Combining (53), (54) and (55) and adding or subtracting terms that are of the 
form (58), (59) or (60) result in the approximation 


1 

eD 

n 




V. - V. 


[1 + J 0 „ ■ 1 ] Q. . + l R. . [N. 
2V X h ij j ij J 



p i ( 



V. 



(61) 
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The only motivation for including a small quantity of the form (60) in the 
above equation is to make the total current (electron plus hole current), as 
predicted by these equations, add up to exactly zero when summed over all 
noninsulated surfaces* 

The important equations for both n and p type materials are listed in 
Table 5-1. 
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Table 5-1. Summary 


denotes the i*-* 1 section of the boundary surface of a quasineutral 
homogeneous region R. The index sets INS, OTH are defined by: 

If ielNS then is an insulated surface 
If ieOTH then A^ is a noninsulated surface 

Other symbols are defined as follows: 


N 

P 


Po 

T 


= Excess electron density 
= Excess hole density 
= Equilibrium electron density 
= Equilibrium hole density 
= Recombination Time 


L 

n 


9 


L 

P 


D 


n» 


D 


P 


V T = KT/e 


e 

K 


= Diffusion lengths for electrons and holes respectively 
= Diffusion constants for electrons and holes respectively 
= Thermal voltage 
= Elementary charge 

= Normal unit vector directed outwards from the region R. 



= Electron, hole, and total currents respectively through the i 
noninsulated surface. These currents are directed 
outwards from the region R. 


£ 


= Position vector representing the point of evaluation of a 
function. 


= Kronecker delta 

Boundary conditions that are assumed for the carrier densities and applied 
potential are: 


If ielNS then the surface A^ has a recombination velocity S^. 

If ieOTH then the surface Aj has constant applied potential and 
constant excess carrier densities denoted by V^, N^, P^. 

All summation symbols denote sums over the index set OTH. All surface 
integrals refer to a unit normal vector (K) directed outwards from the 
region R. 
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Table 5-1. Summary (contd) 


If R is an n Type Material 

For each jeOTH, define the functions 4> j , Yj as solutions to the following 
boundary value problems: 


V 2 <j>. 
J 


0 and V U> . 

J 


h . 

^ in R 
L 

P 


If ielNS and X is on A. then 

i 


^ 6 . • K = 0 and ^ y . 
3 3 


• K = 


-S.x 

l 


2 *3 


If ieOTH and X is on A, then 

l 


<k . = 5 . . and y . = 6 . . . 


To first order in the applied voltages, 


P = p (o) + l V p 

j 


where 


p (0> = I Pj fj and p. = ji- [P. tj - 2 P t V 


Define the constants 


Q. . = J. ^ • ds R. . = X. ^ 4> . • ds 


(Note that Q.. = Q..» R.. = R.. and £ R. . = 0). 
13 x ji 13 ji 13 
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Table 5-1. Summary (contd) 


To first order we have 


p. V. - V. 

= -2 P j [1 + J 2 V t ] Q ij 


(p o + 2 } ^ R ij 

J T J 


Define the constants 


S ijk * 'r *1 ? ? +k ^ T ijk ■ 'r *1 ’ *j • ? ♦k -* 3x 


(Note that T„ k = 
and l S. jk = 0) . 


r [8. . R., + 6.. R. . 
2 ij jk lk ij 


6 jk R ik ] and s kji ■ s ijk + [5 jk - V R ik 


To first order we have 


V. - V. 


-( 2 + n o + D 2 P i ) ? vt R ij + ? (N j “ D 2 P j ) [1 + 2V ^ 1 R ij 

n jT J j J n J T J 


D V. - V. D V. 

+ B 2 J P j [1 + J 2V ] Q ij + 2 D 2 2 V 1 P k S iki 
nj J T J n j,k V T 2 


— = -(^ D N. + D n + f D P. + D p ) 7 ^ R. 
e 2 n i no 2pi p r o r i 


Pi P 0 i xi 1J 
J V *|» 


V, - V, 


+ 7 (D N. - D P.) [1 + — 1 ] R. . + 2 l ~ D P. S... 

4 n j p 3 2V„ xj A. V„, p k lkj. 


2V T « jtk V T P k ife j' 
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Table 5-1. Summary (contd) 


If N^, P^, << n Q , an approximation for the current is 


T 2 2 

V IT V. - V. IT V. 

5F I p j [1 * + i R ii [N j * ?2 f i 

n J J n 


If R is a p Type Material 


For each jeOTH, define the functions 4>j, Yj as solutions to the following 
boundary value problems: 


V 2 * 0 and V 2 «|». = — ^ in R 

2 J L 


If i e INS and X is on then 


V <J>. • K = 0 and v <|r. • K = — i|i. 


2 T j' 


If i e OTH and X is on then 


<b . = 6 . . and f . = S . . 
j ij 3 iJ 


To first order in the applied voltages. 


n = n (o) + l v n, 
j J J 
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I 


Table 5-1. Summary (contd) 


where 


(o) 


I N. and n j = 2Vl [ *j \ N k ♦lc ‘ N j V 

1 IK 


Define the constants 


Q. . = S A $ • ds R. . = J" ^ • ds 

ij A. T J ij A £ T j 


(Note that Q.. = Q . . , R.. = R.. and 7 R. . = 0) . 
ij Ji iJ Ji j iJ 


To first order we have 


n. V. - V, N V 

eiT = 2 N j [1 + 2V X 1 Q ij " (n o + 2“ ) ^ V T R ij * 


Define the constants 


S ijk - 'r *i 5 tj ' ? d3x T ijk - J-R *i « *j • ’ *k ^ 


(Note that : 

and £ S = 0). 
k ijk 


- [6. . R., + 6.. R. . 
2 ij jk lk ij 


6 jk R ik ] and S kji " S ijk + [S jk ^ij 1 R ik 


To first order we have 


eD 


P. D V. D V. - V. 

-(— + p Q + jj- N.) 2 ^ R £j - l (Pj - g- N^ [1 + -Ly— 1 R ij 


D V. - V. D V. 

d" ? N j [1 + X 2 v t J 1 Q ij + 2 5~ A vi N k S ikj 


P 3 


P j »k T 
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Table 5-1. Summary (contd) 


and 


I. , o V 

— = - (£ D P. +D p + r D N. +D 
e 2 p x p *o 2 n 1 n o H V T ij 


V. - V. 
J i 


V. 


- I P| - K N 4 ) [1 + "V* ] R.. + 2 l r 1 D n N k S 


P J n 3 


'H T V T n "k , ifcj' 


If N^, << p Q> an approximation for the current is 


eD 


L V. - V. 

- i I "j [1 * '^2V^ 1 «li 

p J 


L V. 

+IR-4 I- Pi +-f N. (r^- + l) 


ij 


J l 2 "j '2 V t 


V q1 
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SECTION 6 


RECOMMENDATIONS FOR FUTURE WORK 


The analysis presented was intended to be used as an analytical tool for 
the study of latchups but the study was not completed. It is expected that 
this analytical tool will have some useful applications in its present stage 
of development but some important work remains to be done before this tool can 
be applied to the study of latchups. 

Asstiming that the first order perturbation method is valid, the equations 
listed in the summary table have general applicability to homogeneous quasi- 
neutral regions providing the boundary surface is partitioned in such a way 
that p, n and $ are approximately constant on each noninsulated section. What 
is probably the most obvious weakness in the application of these equations is 
the assumption that an entire depletion region boundary surface is a constant 
p, n, <|> surface. The accuracy could be improved by partitioning the deple- 
tion region boundary surface into a few sections with p, n and <|> treated as 
constant on each section. However, such a partitioning would require that we 
properly account for the tangential currents that flow inside the depletion 
region between adjacent sections. It was not shown how to do this and this 
remains to be done. 

A time dependent analysis using low frequency approximations (where 
steady state equations are applied to the interior of the region but DC 
capacitances are associated with various boundary surfaces so that terms 
containing time derivatives of voltages are added to various currents) should 
be relatively straightforward but it was not done here. It is hoped that this 
will be done in future work. 

Examples which illustrate the use of the equations have not been included 
in this publication. It is expected that this will be done in future work. 

Last, but obviously not least, a comparison is needed between measurement 
and theory. This is also expected to be done in future work. 
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Some other assumptions that were made are not essential to the theory 
and are therefore not considered to be a weakness of the theory. These 
include neglecting recombination in the space charge regions, neglecting the 
possibility of avalanching in the depletion regions and assuming that 
junctions are abrupt. These characteristics have no effect on the equations 
listed in the summary table, which is the principal product of this 
publication. These characteristics affect only the relationships between the 
boundary values. Relationships that include such effects as avalanching, 
recombination or graded junctions could be used with the equations in the 
summary table without modification of those equations. 
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APPENDIX A 

HIGH-LOW JUNCTIONS AND DEPLETION REGIONS 


A- 1 


I 


High-low junctions and depletion regions are two examples of the general 
condition of a doping density that is variable in a narrow region of space 
which results in a reasonably well defined space charge region. The analysis 
to follow applies to this general category providing that certain 
approximations are valid. 

We assume that a one dimensional analysis can be used and take the x axis 
to be perpendicular to the junction. Generalized transport equations have 
been provided by Overstraeten, DeMan and Mertens [5]. When applied to the 
interior of the space charge region where the doping density is nonuniform, 
the equations contain an additional drift term and are written as 


J = eD [- T 2 - jp (<(> - DEV)] 

p,x p dx V T dx VT t 


(A-l ) 


J 

n,x 


eD [t 5 - 5- 4- (<t>„ + DEC)] 

n dx V T dx T t 


(A-2) 


where is the total (as opposed to applied) electrostatic potential and DEV 
and DEC are functions of the doping concentrations. The DEV and DEC used here 
differ from those used by Overstraeten et al by a factor of the electronic 
charge. It is assumed that the net (diffusion plus drift) currents for each 
carrier are small compared to the individual contributions [3]. This requires 
negligible recombination in the space charge region. The result is 


d£ 

dx 




- DEV) 


dn 

dx 


n d 
V T dx 


(4> fc + DEC) . 
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A-3 



If the boundaries of the space charge region are at x^ and x 2 » 
can be integrated to give 

P(x,) 

) = exp [ (4> t (x 1 ) - 4> t (x 2 ) + DEV(x 2 ) - DEV^))/!^] 
n(x.) 

n ^ x ) = exp [(<j> t (x 2 ) - 4> t (x 1 ) + DEC(x 2 ) - DEC(x 1 ))/V T ] . 


Let p Q , n Q , <j> Q denote the equilibrium densities and potential, 
and (A-4) to the equilibrium conditions gives 


p (x 2 ) 

~ ^ x ) = exp [(^(Xj^) - 4> 0 ( x 2 ) + °EV(x 2 ) - DEV(x 1 ))/V T ] 


— = exp [ (<|> o (x 2 ) - + DEC(x 2 ) - DEC(x 1 ))/V T ] 


so that (A-3) and (A-4) become 


p(x 2 ) 

pO^) 


p o (x 2 ) 

p ~ (x ) exp “ 4>(x 2 ) )/ v T 3 


n(x 2 ) 

nCx^ 


iTTO ex P [< * <x 2 ) - * (x l» /V T 1 


where 

<j>(x) ■ <|> t (x) - <t> Q (x) 


is the applied potential. To shorten the notation, define B by 


fi = exp [OKx^ - 4 >( x 2 ))/V t ] 


A-4 


the equations 


(A-3) 


(A-4) 

Applying (A-3) 


(A-5) 


(A-6) 


(A-7) 


(A-8) 



so that (A-5) and (A-6) become 


P (x 2 } = P o (x 2 ) 
p(x L ) p 0 (x 1 ) 


(A-9) 


n(x 2 ) _ n 0 (x 2 ) _-l 
n(x. ) " n (x- ) 

1 O 1 


(A-10) 


Let N (x) and N.(x) be the density of donor and acceptor ions. Assuming space 

U A 

charge neutrality at x^ and x^ we get 
p^) - nCxj^) = N A (x t ) - NpCxj^) 
p(x 2 ) - n(x 2 ) = N a (x 2 ) - N d (x 2 ) . 

Combining these equations with (A-9) and (A-10) gives 


P(x t ) 


O' 1 (N D ( Xl ) - n a (x l )1 - [n d (x 2 ) - n a (x 2 )] 

o 1 

p (x ) n (x ) 1 

° 2 B - S' 1 


P» (x l ) 


"o <x l ) 


“ 'W - V x 2» - [ VV - V x i» 

P< * 2> = Pq* 11 !* .-1 n o (x l ) - 

p o <X 2 ) n o <X 2 ) 


(A— 11) 


(A-12) 


A-5 


n(x 1 ) 


[T^f) 13 'W - W 1 " [ W 


W 1 


.-1 . ! o ^4 B 


W 


p o (x l ) 


(A-13) 


n(x 2 ) 


iT5^ B- 1 (n a (x 2 ) - n d (x 2 )J - [N A ( Xl ) - yx^j 


B - S&* .'I 


n o <x 2 ) 




(A-14) 


Since recombination in the space charge region has been neglected, we also have 


J (x. ) = J (x_ ) 
p,x 1 p ,x 2' 


J (x, ) = J (x„) 
n,x 1 n,x 2 


(A-15) 
(A-16 ) 


If the equilibrium carrier densities can be regarded as known, equations 
(A-ll) through (A-16) constitute a complete system of equations describing the 
boundary conditions. But to use these equations in a practical problem, a 
method of estimating the equilibrium densities from existing empirical data is 
needed. Slotboom and De Graaff [4] have related the equilibrium densities to 
the intrinsic densities (the equilibrium densities of the material without 
doping) through the equation 

P o n 0 = n i eXp [AV go /V T ] 

where n^ is the intrinsic carrier density and AV gQ is an empirically determined 
function of the doping level. Empirical values are given by Slotboom and 
De Graaff. Applying this equation to x^ and x 2 yields 


p 0 (x 1 ) n o ( Xl ) = n. exp [AV^x^/V^ 


p o (x 2 ) n o (x 2 ) = n. exp UV go (x 2 )/V T ] 


(A-17) 
(A— 18 ) 


A-6 



If the junction is abrupt with the (known) doping concentration uniform on each 
side of the point of discontinuity, the doping concentrations at x^ and 
are known without the need of solving for x^ and x^ (for a more general doping 
profile, x^ and x^ would have to be solved and their values would depend on 
the doping profile) and AV q (x^), AV q (x£) can be innnediately evaluated from 
the data provided by Slotboom and De Graaff. By assuming that the majority 
carrier density is approximately equal to the doping density, (A-17) and 
(A-18) can be used to solve for the minority carrier densities at x^ and x^. 

The equilibrium minority carrier densities can now be regarded as known. 
Approximating the majority densities with the doping densities may be adequate 
for the purpose of using (A-17) and (A-18) to solve for the minority densities, 
but it is not sufficiently accurate for other calculations for the following 
reason. If the equilibrium majority carrier density is set equal to the 
doping density, the error in the equilibrium majority carrier density will be 
equal to the equilibrium minority carrier density. If this error ultimately 
showed up in the nonequilibrium majority carrier density in the form of an 
additive constant, it would not affect the calculation of the diffusion 
current and it would negligibly affect the calculation of the drift current, 
so the error would be minor. But the error will show up in the nonequilibrium 
density in a way that is more serious than the error resulting from an 
additive constant. This is because an error in the equilibrium density will 
affect the boundary values of the nonequilibrium excess carrier density as 
calculated by (A-ll) through (A-14). This error is small compared to the 
majority carrier density, but it may not be small compared to the excess 
majority carrier density. In fact, it is not difficult to show that the error 
in the calculated boundary value for the nonequilibrium excess majority 
carrier density (as calculated by (A-ll) through (A-14)) exceeds the error in 
the equilibrium density. When the diffusion equations are used to calculate 
the excess density as a function of location, errors in the boundary values 
affect the solution in a nontrivial way and the calculated currents will have 
serious errors. Therefore, to calculate the equilibrium majority carrier 
density, we must be more accurate than simply setting it equal to the doping 
density. A more accurate treatment uses space charge neutrality to get 


( A-20 ) 


P 0 (X 2 ) " n o (x 2 ) = N A (X 2 ) " W * 

Note that this space charge neutrality condition is needed in order for (A-ll) 
through (A-14) to predict the correct equilibrium (13=1) values for the carrier 
densities. 


Example 1: High-Low Junction 

Consider an n - n + junction and let x 2 be on the heavily doped side. 

To use more suggestive notation, let 

*h - V “ft - *1* N Dh “ W' N Dft • N D (x 1 ) 

p oh = p o (x 2 )> p oft * p o (x l )> "oh * "o (x 2 ' ’ "oft * "o <x l ) 
p h = P(x 2 ), P ft = pCxj^), n h = n(x 2 ), n ft = n(x^ 

4> h ■ ♦(* 2 )» ♦jt - ♦(*]_) • 

Equations (A-ll) through (A-14) become 

(A-21 ) 


(A-22) 


(A-23) 


— B N Q - N 

_ n 0 ft D11 Dh 

p 2. “ p . n . . 

*oh _ oh „-l 


p 0 a 


n 


oi 


Ph - 


_oi 

'oh 


S N Dh " N DH 


— B" 1 


ofi. 


■oh 


n 


oh 


n Sl " 


^oh 

PqI 


B N DH + N Dh 


oh „-l 


fi - 


oi 


■oh 

? oi 


l 
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n h - 


— B 1 N + N 0 
P oh Dh Dft 


n 


(A- 24) 


oft 


n 


B - 


■oft „-l 


oh 


oh 


where 


B = exp - <|> h )/V T ] . 


(A-25 ) 


Assume that AV^ q is essentially zero on the low doped side and represent its 
value as simply AV^ q on the high doped side. Equations (A-17) and (A-18) 
become 


P oi N Dfe ' ”l 


P oh N Dh * "i • 


Combining these equations with (A-19) and (A-20) yields 


2 

v - n db + 4 


(A-26) 


2 

n oh = N Dh + eXp tAV go /V T ] ‘ 


(A-27 ) 


Using the same analysis for a p - P + junction yields 


P ft = 


r: B_1 n ai + % 

oft 


■oh 


n 


oh „-l 


B - — B 
P oft n oft 


A-9 



n 


oil 


n 


p h = 


oh 


6 N Ah + 


- 


n 


oft 


■oh 


n 


oh 


n„ = 


■oh 

'oft 


13 N Aft " N Ah 


n 


oh „-l 


n 


B 


oft 


■oh 

'oft 


n h = 


— fl 1 N - N 0 
P oh Ah Aft 


n 


oft 


oh 


B - 


■oft n -l 


•oh 


n 


2 

n. 

l 


oft N 


Aft 


2 

n. 

n = rp- exp [AV /V ] 
oh N Ah go T 


2 

n. 

P oft = N Aft + n7. 

Aft 


2 

Poh ■ N Ah + ex P W V go /V T ] 


Example 2: PN Junction 

Let x^ be on the p side and x^ be on the n side. To use more suggestive 
notation, let 
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X p - V X n ■ V N Ap ’ “a* 31 ! 1 - N Dn = N D (x 2 ) 


P or1 = P (x..), ? n „ = P (x ), n = n (x ), n = n(x) 
op o 1 on o l op o 1 on o 2 . 


op 


P p = P^), P n = p(x 2 ), n p = n(x 1 ), n n = n(x £ ) 


+p = + (x l ) * *n = <{,(X 2 ) ’ AV p = ^go^* AV n = AV go (x 2 } 


Equations (A-ll) through (A-14) become 


P„ = 


n i 

on — 1 __ 

13 N a + N 
n Ap Dn 

_22 

n i P 

_on B -1 . Ion B 

n op p op 


p n = 


B N. + N. 
n Dn Ap 

on 

P i n 

-ae b- 1 - -52 s 

n 


on 


on 


n = 


— B N + N 
p AP Dn 

122 


P n . 


on „-l 


n 


B 


op 


on 


op 


8 


n = 
n 


^ B 1 N_+ N. 
P on Dn Ap 


P i n 

122 d -1 _ _22 


on 


n 


on 


with 


B = exp [(<l> p - <l> n )/V T ] 


A-ll 


while (A-17) through (A-20) imply 


2 

p on ■ 5T exp [4 W 

Dn 


2 

n. 

n = rp exp [AV/V] 
op N Ap p T 


2 

n. 

p nn = N a + ST" ex P UV/Vj 
op Ap N a p p T 


2 

"on = N Dn + 5T ® Xp [A W * 
Dn 
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APPENDIX B 

A UNIQUENESS THEOREM 


B-l 


Let a region R be bounded by a surface A. Consider the boundary value 
problem 

(V 2 - K)f(X) = f(X) for X e R 

a(X)«|»(X) + b($) $Y(X) • n = g (i) for X e A 

where f and g are specified functions and n is the outer normal unit vector. 

K is a constant satisfying 

K > 0 . (B— 1 ) 

At any point X on the surface A, the functions a and b are to satisfy one of 
the following conditions: 

Either 


a(X) = 0 and b(X) * 0 


(B-2) 


or 


a(X) 0 and b(X) = 0 


(B-3) 


or 


a(X) 4 0 and b(X) t 0 and adt)/b(X) > 0 . (B-4) 

It will be shown that the solution f is unique* Note that this is not the 
most general possible uniqueness theorem because of the conditions (B-2), 

(B-3), and (B-4) (there are also uniqueness theorems, for example, to the 
Helmholtz equation, which does not satisfy (B-l)) but it is general enough for 
our purposes since recombination velocities and diffusion lengths are positive 
quantities . 
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B-3 


Let »|»^ and ^ be two solutions to the boundary value problem and let 

= ^2 " ♦]_ • (B-5) 

Note 


(V 2 - K) f = 0 for X e R 


a(X) <|»’(X) + b(X) ^ H»'(X) • n = 0 for X e A 


(B-6) 
( B— 7 ) 


Green's first identity states that 


/ u,' V 2 y d 3 X + /„ $ Y’ • $ y d 3 X - J*. • ds . 


2 

Using (B-6) to replace V i|/' in the first integral gives 


K X R (»|»') 2 d 3 x + X R |$ y I 2 d 3 X = J A f ^ • ds 


(B-8) 


Let A q denote the section of the surface A such that either (B-2) or (B-3) is 
satisfied. Note that (B-7) combined with either of these conditions implies 
that 


Y' V • n = 0 for X e A . 


Therefore, 


X^ f ^ y' • ds = / A , <J»' $ <J»’ * ds 


where A' is the section of A such that (B-4) is satisfied. Using (B-7) gives 


/ f' $ t* • ds - -X A , ^ W ) 2 


ds 


b(X) 


B-4 



so that (B-8) becomes 


K J R (f) 2 d 3 X + J R |^'| 2 d 3 X + X A , g (f) 2 ds = 0 . 

From (B— 1) and (B— 4) it is seen that the left side of the above equation is a 
sum of nonnegative quantities. Therefore, the equation implies 

Y' = 0 

or, from (B-5), 

*1 - *2 • 


B-5 



